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Abstract 

We consider systems governed by the Gross Pitaevskii equation (GPE) with the Morse potential V(x) = —2e“““) 

as the trapping potential. For positive values of the coupling constant g of the cubic term in GPE, we find that the 
critical value gc beyond which there are no bound states scales as (for large D). Studying the quantum evolution 
of wave packets, we observe that for g < gc, the initial wave packet needs a critical momentum for the packet to escape 
from the potential. For g > gc, on the otherhand, all initial wave packets escape from the potential and the dynamics is 
like that of a quantum free particle. For g < 0, we find that there can be initial conditions for which the escaping wave 
packet can propagate with very little change in width i,e., it remains almost shape invariant. 
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1. Introduction 

The wave packet dynamics of the usual Schrodinger 
equation has been extensively studied for the free parti¬ 
cle as well as for different kinds of external potential. For 
the free particle, if a nonlinear term in the wave function 
is added to the Schrodinger equation, then one has the 
non-linear Schrodinger equation (NLSE)[i|, [^. This too 
has been a subject of intense study and one has found 
various kinds of exact solutions like solitons, pulses, fronts 
etc.[i,|3,[i|. If one adds a potential, then one has the Gross 
Pitaevskii equation (GPEjig which has been very use¬ 
ful for describing the process of Bose Einstein condensate 
(BEC)[1, 0 for a trapped gas (it should be noted that 
experimentally BEG has only been observed in a trapped 
gas). The trapping potential makes the problem far more 
difficult and the only case which has been reasonably well 
studied is the simple harmonic potential. In an experi¬ 
ment, it is actually far more convenient to tune the trap¬ 
ping potential and observe different behaviors of the con¬ 
densate. With this in mind, we wanted to study the GPE 
with a Morse potential so that one can achieve a much 
greater flexibility in adjusting the potential. In fact, the 
potential can have bound states or no bound state at all, 
depending upon certain parameters of the potential. The 
corresponding dynamics is treated both for g > 0 and 
g < 0, where g is the coupling constant of the non linear 
term. 
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The primary importance of GPE lies in the fact that 
it describes the dynamics of the condensate in the pro¬ 
cess of Bose Einstein Condensation (BEG) [13, [HI. The 
emergence of different kinds of soliton (dark, bright, grey 
12|- 1^ is well known in case of GPE. Also the control¬ 
lable soliton emission [l^ has been investigated for GPE 
with a shallow trap and with negative interspecies interac¬ 
tion. In this article, we have dealt with the attractive and 
repulsive interatomic interaction separately for the Morse 
potential and the dynamics in both of these two cases have 
been explored both numerically and analytically. 

Time dependent Gross Pitaevskii equation (GPE) with 
one spatial dimension has the following form: 


ihdt’iji = + Vext{x)'tp 

2m 


( 1 ) 


2. Bound state 

The stationary states of Gross Pitaevskii equation (GPE) 
have the form'0(a;,t) = sothat^u(a;) = — 

gu\u\'^ -\- Vextu and lowest value of p is the ground state 
energy Eq which is to be obtained by minimizing E of Eq. 

@. 

/ (x> j-2 „ 

[—\\7ip\^ + Vext\ili\^ + ^\ilj\^]dx ( 2 ) 

where, Vext{x) = D{e-^^^-2e-^^). For o « 1, Vext{x) - 
x^ indicating the oscillator limit of the Morse. On the oth¬ 
erhand, for a -A oo, Vext{0) = —D and 0 else where for 
a; > 0 as is clearly depicted in Fig ©• 
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Figure 1: Schematic diagram of Morse potential depicting the two 
extreme conditions. Blue curve shows approximately the oscillator 
limit (a <C 1) and the violet one (for a oo) indicates the free 
behavior for region of a: > 0. 


With the differential equation for u given as 

we can explore the asymptotic solution as a; —> — oo. The 
function u has to vanish as x —> ±oo being a bound state 
function and hence the dominant part of Eq. for X —>■ 
—oo will be and we have in that range 


(Pu ^ ^ 
2m dx^ 


—2ax 


U 
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( 4 ) 


Making the substitution = y, we find the exact asymp¬ 
totic form of u for X —>■ — 00 (y —>■ oo) to be of the following 
form 


Uasy = Ae (5) 

where, and A is a constant to be determined 

from normalization in y-space. Since the ground state 
wave function will have no zeroes at any finite nonzero 
value of y (note that — oo < x < oo corresponds to 0 < 
y < oo), we take u{y) = y/Ny°‘e~^y, to be the trial wave 
function for the minimization of E given in Eq. m and 
we obtain the following expression for the ground state 
energy, 


E{a) 



■^a{l-2K) + 


agm r(4a) 
fi2a2 24“[r(2a)]2J ^ ’ 


If g = 0 (i.e., usual Schrddinger equation with a Morse 
potential), we get the critical a for minimization of energy 
is a = AT — ^ and the ground state energy turns out to be 


which is the exact answer.lt should be noted that for K = 
the binding energy reaches zero, i.e., for A' < 4 (equiv¬ 
alently for a > -\J , there can be no bound state in this 

potential. We now analyze the effect of non linear term. 
Returning to Eq. ®, we consider the term involving ‘gi’ 


and and writing A = ff^, we get as the minimization 
condition 


2a-f (I - 2A:)-f A 


d r(4a) 

Sr 24“[r(2a)]2 


= 0 


( 8 ) 


Since the tightly bound situation correspond to large value 
of AT, we use the above equation in the limit when a is rea¬ 
sonably large and r(x) can be replaced by the asymptotic 
form 

r(x) ^ e~^x^x~^^'^\f2^\V -|- —I-] 

We get from Eq. (|5]) to the leading order, a to be given 

by 

y. 1 ^ 1 

“ 2 2y^(2if-1)1/2 

This gives a ground state energy of 




an equation which is accurate only for AT ^ 1. In Table 1, 
we show the comparison of ground state energies obtained 
from the actual solution of the minimization condition of 
Eq. ([5]) and the energy obtained from Eq. ® for fixed 
values of A of A = 1.0 < Ac and different values of A'. As 
expected the approximation of Eq. © gets closer to the 
answer obtained from Eq. ® as AT increases. The increase 
of A weakens the ground state and for large AT, we have 
the result that the critical A for disappearance of a bound 
state scales as AT^/^ i.e., H^/^. 


K 

Eo{fromEq.^) 

Eo{fromEq.^) 

2 

-1.762 

-2.243 

3 

-5.619 

-6.25 

4 

-11.504 

-12.25 

5 

-19.404 

-20.25 

6 

-29.315 

-30.25 


Table 1: Results compairing the ground state energies from varia¬ 
tional calculation (third column) with the expression (Eq. l(9} ob¬ 
tained considering the asymptotic series for different values of K 
keeping A fixed at 1.0 


The interesting thing would be to explore the dynamics 
of an initial Gaussian wave packet in the region A > Ac and 
A < Ac for A > 0. In the next section we describethe wave 
packet dynamics both for A > 0 and A < 0. 


3. Wave packet dynamics 

The effective Hamiltonian of GPE with a Morse trap 
is given by 

H=^+ - 2e-“") +gj dx\^{x)\'^ 

We make E[ dimensionless by the following rescallings: 
p = ^ X = f, A = f, t = 


2 
























V = ^ and —= 7 —. Hence the relation between A 

and 7 turns out to be A = Considering these 

transformations the dimensionless Hamiltonian turns out 
to be 




2 e “) + •\/^7 


dx\'ijj{x)\‘^ (10) 


We consider Gaussian wave packets of the form given be¬ 
low 




1 


(»-xo(t))^ 

2A(t)^ 


^ip{t)xlh 


( 11 ) 


i.e., we assume an initially Gaussian form remains Gaus¬ 
sian with a changing centre and width. The initial shape 
corresponding to So = 0 and A = Aq has the energy 


(Eq) 
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1 

2^2 Ao 



[e' 


A" 



( 12 ) 


The dynamics of the wave packet is governed by the fol¬ 
lowing equations (we drop the bars with the understanding 
that all quantities are dimensionless). 


d{x) 

dt 

d{x'^) 

dt 

d{p) 

dt 

d{p^) 


= (P) 

= {xp + px) 

= -(^) = 2 ((e-^^) - (e-^)) 
dx 


, dV dV , /— d 

= -{P^ + -rP/ - ^27r7— 
dt dx dx dt 


d{xp + px) 
dt 


= 2(p2) - 2aD{x 


dx 


\ij}\^dx (^^) 

) -I- J IV'I"' 


dx 


= 2(p^) -I- 4aiA((a::e 
^dx 




xe-n) 


+ 


V^7 J IV'I"' 


With the help of the equations given in Eq. (USD, we find 
that the dynamics of the wave packet will be governed by 
the two coupled equations given below 

^xo = (14) 



2 1 


-f -h 4A2 

A 



_ g-(2a:o-A^)j 


(15) 


Here we have explored the wave packet dynamics for K = 2 
which fixes 7 c = 0.917. 


3.1. Dynamics: (7 > 0 and 7 < 7 c); Corresponding 
figure: Fig[^ 
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Figure 2: Figure (a) shows the dynamics of the peak (fco) of the wave 
packet for 7 (= 0.5) < 7c(= 0.917) for four different values ofpo- The 
potential parameter K = 2 and the initial xq = 0 while the width 
Ao = 0.4. For low values of po the wave packet hits the potential 
barrier at right and reflects back and forth within the trap showing 
oscillatory behavior. With increase of po the number of reflection 
reduces and finally at po = Pth the wave packet comes out of the 
potential barrier which is represented by the linear increase of xo{t) 
at later time (green curve). The figure above shows that in this case 
p^^ = 0.45. Figure (b) shows the dynamics of the width (A). All the 
parameters are same as are taken in figure (a). With lower values 
of Po when the wave packet hits the potential barrier and reflects 
back and forth within the trap, A{t) also shows nearly oscillatory 
behavior. With increase of po the the wave packet finally escapes 
with increase of A(t) as t increases (green curve). 

In this parameter region after solving the coupled equa¬ 
tions given in Eq. (1141) and Eq. m numerically, we have 
found that when initial momentum (po) of the wave packet 
is small, it reflects back from the potential barrier at right. 
Then it moves to the left and collides with the infinite bar¬ 
rier or the potential at right. This back and forth oscil¬ 
lation within the trap continues with time. With the in¬ 
crease of Po, the oscillation of the peak of the wave packet 
within the trap goes on with the higher amplitude up to a 
certain threshold value {{pth) of po. At po = pth it simply 
comes out of the potential barrier and hence the oscillatory 
behavior of a:o(t) ceases. Eor po > Pt) the graph for xo(t) 
vs. t shows sharp linear increase (not shown in given fig¬ 
ures). Eor K = 2, we observed pth = 0.45. Hence with the 
help of Eq. (HD), we conclude that the initial wave packet 
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has to have the minimum average energy E = 0.751Z? to 
come out of the potential barrier. We will call this mini¬ 
mum average energy required for emitting the wave packet 
from the potential as the threshold energy {Eth)- It should 
be noted that in the classical case, for a partice to es¬ 
cape from the Morse potential, the average E needs to be 
greater than or equal to zero. If the particle is at the origin 
initially with a potential energy of -1, then the threshhold 
momentum would be po = \/2. In the quantum case the 
average initial momentum is 0.45 clearly showing the role 
of quantum fluctuations. The dynamics of the width also 
follows the same qualitative feature as described above. 
The dynamics in this parameter region is clearly depicted 
in Fig [21 For convenience we keep the initial width of the 
wave packet fixed at Aq = 0.4 throughout the study. 

3.2. Dynamics: (7 > 0 and 7 > 7 c); Corresponding 
figure: Fig [3] 



Figure 3: Figure (a) shows the dynamics of the peak (xq) of the wave 
packet for 7 = 1.2 for three different values of po- The dynamics of 
xq always shows linear increase for all these three values of po unlike 
the previous case of 7 < 7 c. This implies in this parameter region the 
wave packet will always come out of the potential barrier whatever 
be the initial momentum. The width also increases and the dynamics 
is like that of a free particle. 

In this parameter region we have observed the linear 
increase of Xo(t) irrespective of initial momentum (po) of 
the wave packet. The width also shows continuous increase 
with time. With the help of Eq. m, we now conclude 


that wave packet in this parameter region always have en¬ 
ergy E > Eth{= 0.7511? obtained ioi K = 2 and 7 = 0.5). 
Thus the wave packet always comes out of the potential 
barrier and delocalises in space. In Fig [31 we clearly de¬ 
scribe this feature considering j{= 1.2) > 7 c(= 0.917), 
Aq = 0.4 and K = 2. The behavior is always thatof a free 
particle. 

3.3. Dynamics (7 < 0); Corresponding figure: Figj^ 
and Fi^ 



Figure 4: Figure (a) shows the dynamics of the peak (xq) of the 
wave packet for 7 = —0.5 for three different values of po- K = 2 and 
(Aq = 0.4) are considered. Like Figure 2(a), with lower values of 
po the wave packet hits the potential barrier and reflects back and 
forth within the trap. xo(t) shows oscillatory behavior in this case 
also. With increase of po the number of reflection reduces, oscillatory 
behavior ceases and finally when po = Pthi= 1-35) the wave packet 
comes out of the potential barrier indicated by the gradual increase of 
xo(t). In figure (b), the dynamics of width shows a random oscillation 
for po < Pthi= 1-35) and as po —Ptht oscillation takes more 
periodic manner with constant amplitude. But unlike the previous 
cases here the width is bounded by a maximum and minimum value 

(Ao). 


In this parameter region the dynamics of xq shows the 
same qualitative feature as is described in Figure [2] We 
set 7 = —0.5, K = 2 and Aq = 0.4. But unlike the 
previous cases now the dynamics of the width is bounded. 
The random oscillation of width continues with time with 
moderate amplitude and as po —>■ pth the oscillation takes 
place in a more periodic manner with constant amplitude 
which implies that the wave packet will never delocalise 
in space after coming out of the potential barrier. This 
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case is depicted in Fig [4] In Fig [5] we have shown the 
dynamics of xq and A for a greater value of 7 = —1.2. 
The dynamics of xq follows the same qualitative feature 
as described earlier. But this time the dynamics of the 
width shows an interesting feature. The width of the wave 
packet has a tendency to retain its initial shape within the 
trap as well as after coming out of the trap. 



Figure 5: Figure (a) shows the dynamics of the peak (xq) of the 
wave packet for 7 = —1.2 for three different values of po- K = 2 
and (Ao = 0.4) are considered. Like Fig. 4(a), with low values of 
po the wave packet hits the potential barrier and reflects back and 
forth within the trap, showing oscillatory behavior. With increase 
of po the number of reflection reduces, oscillatory behavior ceases 
and Anally when po = Pt/i(= 1.35) the wave packet comes out of 
the potential barrier indicated by the gradual increase of a:o(t)- In 
figure (b), the dynamics of width shows periodic oscillation for all 
values of po and as po Pth-, fhe amplitude of oscillation reduces. 
At Po = Pth-> the width remains almost constant for all time. 

In Table 2., we have calculated and compared the values 
of threshold energy for different values of K. For 7 > 0 , 
as K increases the value of Ett increases which is justified 
by the fact that with increase of itT, Morse potential ap¬ 
proaches oscillator limit (more highly trapped limit) and 
consequently greater energy is required for escaping from 
the potential. 


K 

Ethiior 7 = 0.5) 

Ethiior 7 = -0.5) 

2 

0.75ID 

-0.934D 

3 

0.904D 

-1.09D 

4 

0.97ID 

-1.I4D 

5 

0.986D 

-1.I7D 

6 

1.004D 

-1.I85D 


Table 2: Results showing how the threshold energy (Etji) behaves 
with potential paramater K for two fixed values of 7 (7 = 0.5 and 
7 = -0.5) 

4. Conclusion 

In conclusion, we have explored the features of the GPE 
with the Morse potential. We have found that for a posi¬ 
tive coupling constant and a deep potential there is a crit¬ 
ical coupling gc where the ground state disappears as a 
bound state and scales as the depth of the potential. 
For the dynamics in this potential if 5 < gc-, the initial 
wave packet needs to posses a threshold average momen¬ 
tum for the packet to escape from the potential. For g > gc 
however, the wave packet dynamics always resembles that 
for a quantum free particle. If the coupling g is negative, 
we find that the packet does escape from the potential for 
above threshold value of the average momentum and fur¬ 
ther if g is more negative than a critical value, the width 
of the packet remains constant in time. 


Acknowledgments 

One of the authors, Sukla Pal would like to thank S. N. 
Bose National Centre for Basic Sciences for the financial 
support during the work. Sukla Pal acknowledges Harish- 
Chandra Research Institute for hospitality and support 
during visit. 

References 

[1] N. Gisin and M. Rigo, J. Phys. A: Math. Gen. 28 (1995) 
7375-7390. 

[21 A. S. Fokas and A. R. Its, J. Phys. A: Math. Gen. 37 (2004) 
60916114. 

[3] Vivek M Vyas, J. Phys. A: Math. Gen. 39 (2006) 9151-9159. 

[4] W. van Saarloos and P. G. Hohenberg, Physica D, 56 (1992) 
303-367. 

[5] W. van Saarloos and P. C. Hohenberg, Phys Rev Lett, 64 (1990) 
749-752. 

[ 6 ] E. P. Gross, Structure of a quantized vortex in boson systems, 
Nuovo. Cimento., 20 (1961) 454-457. 

[7] L. P. Pitaevskii, Soviet Phys. JETP, 13 (1961) 451-454. 

[ 8 ] K. B. Davis et al. Phys Rev Lett 75 (1995) 3969. 

[9] A. Leggett Rev Mod Phys 73 (2001) 307-356. 


5 











[10] Franco Dalfovo, Stefano Giorgini, Lev P. Pitaevskii, and Sandro 
Stringari, Rev. Mod. Phys. 71 (1999) 463. 

[11] C. J. Pethick and H. Smith. Bose-Einstein Condensation in 
Dilute Gases. Cambridge University Press, Cambridge, 2002. 

[12] S. Burger et al. Phys. Rev. Lett, 83 (1999) 5198. 

[13] P. G. Kevrekidis et al. Phys. rev. Lett, 110 (2013) 264101. 

[14] Th. Busch and J. R. Anglin, Phys. Rev. Lett. 84 (2000) 2298. 

[15] Mara 1. Rodas-Verde et al. Phys. Rev. Lett. 95 (2005) 153903. 


6 



